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Singular Charge Density at the Center of the Pion? 
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We relate the three-dimensional infinite momentum frame spatial charge density of the pion to its 
electromagnetic form factor F 7r {Q 2 ). Diverse treatments of the measured form factor data including: 
phenomenological fits, non-relativistic quark models, the application of perturbative QCD, QCD 
sum rules, holographic QCD and the Nambu Jona-Lasinio (NJL) model all lead to the result that 
the charge density at the center of the pion has a logarithmic divergence. Relativistic constituent 
quark models do not display this singularity. Future measurements planned for larger values of Q 2 
may determine whether or not a singularity actually occurs. 
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Understanding the pion is a necessary step to learning how QCD describes the interaction and existence of ele- 

1 ' mentary particles. As a nearly massless excitation of the QCD vacuum with pseudoscalar quantum numbers, the 
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pion plays a central role in particle and nuclear physics as a harbinger of spontaneous symmetry breaking and as the 
carrier of the longest range force between nucleons. 
— i , The importance of the pion has been recognized by a huge level of both experimental and theoretical activity aimed 
at measuring its properties and understanding its structure. New measurements of the pion electromagnetic form 
factor, F 7T (Q 2 ), have been performed [J0 and are planned Q. Here we present the first phenomenological analysis 
q ■ of existing data to determine the charge density of the pion in a model independent manner. 

A proper determination of a charge density requires the measurement of a density operator. We shall show that 
measurements of the pion form factor directly involve the three-dimensional charge density of partons, in the infinite 
momentum frame, poo(a;~,b). In this frame [J], the electromagnetic charge density J° becomes J + and 

Poa(x~,b) = J+0~,b) = ^2e q q(x~ ,b)-/ + q(x- ,b) = e q V2q\(x~ , b)q+(x~ , b), (1) 
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where q + (x fJl ) = 7°7 + /v / 2g(a; A1 ), the independent part of the quark-field operator q{x IJ '). We set the time variable, 
x + = (t + z)/y2, to zero, and do not display it in any function. 

We are concerned with the relationship between charge density and the electromagnetic form factor F^^Q 2 ), which 
is determined from the current density via the relation: 

2 (y + , P y+(o)|p+,p) 
KKQ) = ^+ > V) 
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where states are normalized as (p' + , p'|p + , p, ) = 2p + (2ir) 3 6(p' + — p + )6^ (p' — p). We take the momentum transfer 
la = p' a — Pa to be space-like, with Q 2 = —q 2 > 0, and use the Drell-Yan (DY) frame with (q + = 0, Q 2 = q 2 ). The 
matrix element appearing in Eq. Q involves the combination of creation and destruction operators: 6 T 6 — Sd for 
each flavor of quark, so that the valence charge density is probed. Note also that the form factor F\ is independent 
of renormalization scale because the vector current q^^q is conservedjil. 

The spatial structure of a hadron can be examined if one uses H @ 0] states that are transversely localized. The 
state with transverse center of mass R set to is formed by taking a linear superposition of states of transverse 
momentum: 



|p+,R = 0,A)^AA J -fl>\p+,p,\), 



(3) 



where |p + ,p,A) are light-cone helicity eigenstates j§} and N is a normalization factor satisfying |A/"| 2 / wj^jz = 1- 
Wave packet representations can be used to avoid states normalized to S functions [l(| , but this leads to the same 
results as using Eq. ([3]). Considering that 2p + p~ — p 2 = m 2 > 0, one finds that p + must approach infinity. This 
ultra-large value of p + (infinite momentum frame) maintains the interpretation of a pion moving with well-defined 
longitudinal momentum[9(. It is in just such a frame that the interpretation of a hadron as a set of a large number 
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of partons is valid. Setting the transverse center of momentum of a state of total very large momentum p + to zero as 
in Eq. ©, allows the transverse distance b relative to R to be defined. 
Next we relate the charge density 

, _ . , _ (p + , R = 0, Aj poofr-, b) \ P +, R = 0, A) 
Pool2: ' j ~ (p+,R = 0,A|p+,R = 0,A> ' { ) 

to F n (Q 2 ). In the DY frame no momentum is transferred in the plus-direction, so that information regarding the x~ 
dependence of the distribution is not accessible. Therefore we integrate over x~ , using the relationship 



to find 



ql(x-,b)q+(x-,b) =e i ^ x ~e- i ^ b ql(0)q + (0)e i ^ b e- i ^~, (5) 

P (b)= J dx- /9oo (a;-,b) = (p + ,R = 0,A|p oo (0,b)|p+,R = 0,A) /(2p+). (6) 

Furthermore, the use of Eqs. (15I3I2[) leads to the simplification of the right-hand-side of the above equation: 

P(b) = J ^^(Q 2 = q 2 )e-^ b , (7) 

where p(b) is termed the transverse charge density, giving the charge density at a transverse position b, irrespective 
of the value of the longitudinal position or momentum. This relation between an integral of the three-dimensional 
infinite momentum frame density and the electromagnetic form factor is our principal new formula. Previous results 
0, H S Eli 03 involved the integral over the longitudinal momentum fraction x of the impact parameter parton 
distribution function (pdf) q(x, b), which gives the charge density for a quark at position b for a momentum fraction 
(of the plus-component) x. The equality of the respective integrals over x~ or x of the quantities poo(x~ , b) and q(x, b) 
is an example of Parseval's theorem. The central charge density of the pion is determined by p(b — 0), because the 
longitudinal dimension is Lorentz contracted to essentially zero in the infinite momentum frame 

Recent pion data[l|, l2| provide an accurate measurement of the pion form factor up to a value of Q 2 — 2 Ah GeV 2 . 
Their analysis includes an assessment of the influence of the necessary model dependence caused by extracting the 
form factor from the measured cross sections on the experimental error bars. The existing data for the pion form 
factor show that it is well represented by the monopole form 

FAQ 2 ) -V(i + R 2 Q 2 /6), (8) 

with R 2 = 0.431 fm 2 . A better representation of the data may be a monopole plus dipole 0] which involves the square 
of the term of Eq. ((8]) , but any form involving the monopole term leads to a singular central charge density. This is 
because the use Eq. ([5]) in Eq. ([7]) gives the result: 

P( b ) = 75^' 9 

where Kq is modified Bessel function of rank zero. For small values of b this function diverges as ~ log(6). This 
divergence is very surprising because the charge density we are considering measures a valence quark operator between 
eigenstates of the full Hamiltonian. The divergences of quark distribution functions that occur at small values of 
Bjorken x do not occur here. Any model, such as vector meson dominance or holographic QCD [l3|, [l4|, HH that yields 
a monopole form factor has a central density with a logarithmic divergence.. 

Intuition regarding a possible singularity in the central charge density may be improved by considering other 
examples. Suppose that the non-relativistic (NR) limit in which the quark masses are heavy is applicable. In this 
case, the pion would be a pure qq object and the charge density is the Fourier transform of the form factor. Given 
the form factor of Eq. ||5J) the three-dimensional density is uniquely given by 

. . 3 -VSr 

pNR[r) = - « 10 

2 7r r K 



where r is the distance relative to the pion center of mass. If one takes r = \Jb 2 + z 2 as demanded by the rotational 
invariance of the non-relativistic wave function, then one finds dzpNn(r) is equal to p{b) of Eq. ([5]). This is 

expected because in the NR limit the charge density is the same in all frames, including the infinite momentum 
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frame. The meaning of the X/r behavior of the density can be understood by considering that for a qq pion, the 
wave function is the square root of the density so that the short distance wave function j/j^r ~ 1/\A"- Using the 
Schrodinger equation, one finds that the potential must contain terms proportional to 1/r 2 . There is no evidence that 
the strong interaction potential behaves in this manner. Thus a non-relativistic viewpoint tells us that the central 
singularity derived from the form factor falling as 1/Q 2 requires unsupported assumptions regarding the nature of 
the short-distance interactions between quarks. On the other hand, the lowest-energy solution of the Dirac equation 
for hydrogenic atoms has a singular radial behavior, ^£>(r) ~ l/r 1- ' 1 'e -r (7 = vl — Z 2 a 2 ), near the origin at r = 0. 
Consider puip) = dz|^£)(r)| 2 and define 77 = 2 — 27, which ranges between and 2. We find for small values of b 
(in units of twice the appropriate Bohr radius) that po{b) is well behaved for < 77 < 1, behaves (for all b) as K Q (b) 
for 77 = 1, and behaves as l/W for 1 < 77 < 2. Thus there are physical examples with a singular central density. 

The divergence of the central transverse charge density encountered here may be the consequence of using a simple 
parametrization, so we shall consider the predictions of a variety of different approaches. We begin with perturbative 
QCD (pQCD) which provides a prediction [l6|, [l7| for asymptotically large values of Q 2 that 

Um Q 2^ 0O F 7r (Q 2 ) = 16™ s (Q 2 )/ 2 /Q 2 , (11) 

with the pion decay constant f n = 93 MeV, and in leading order: 

47T 

(ll-|n/)In^ 

with n f the number of quarks of mass smaller than Q and A is a parameter fixed by data. One might think that the 
logQ 2 term in the denominator would lead to a non-singular behavior of p(b) for small values of b. This is not the 
case. To see this, consider the integral: 

dQ MQb), (13) 



QlogQ/A 

for the case Q ma x > Qo > A > 0. In the limit that (5 max approaches infinity, this is the contribution of the integral of 
Eq. arising from values of Q > Q , assuming that the value of Qo is large enough for Eq. IT2|) to be valid. Take 
Qmax = where e is a small positive number such that Jo(e) = 1 to any specified degree of numerical precision. 
Then 

C qt^m Jo(Q6) = loglog( ^ ) - loslos( x ) ' Q - b « 1 (14) 

In the limit that b approaches zero Eq. (fT4|) becomes 
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We see that the pQCD form factor corresponds to a singularity at short distance. The same feature would arise in 
any model form factor such as those based on sum rules e.g (l8j that joins smoothly to the pQCD result at very large 
values of Q 2 . 

Chiral quark models (see the review [lj|) present other examples of transverse charge densities that are singular at 
the center. In those models, the pion form factor takes the monopole form of Eq. ([8]) so that the central density diverges 
as log b at the origin. Nevertheless all physical observables, including f„ and structure functions, are computed to be 
finite. We consider two such models. The first is the spectral quark model SQM [2(1 HH In this model F n (Q 2 ) takes 
the form of Eq. ((5]) with i? 2 /6 = m 2 The impact parameter dependent parton distribution function is given by (2lj 
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q(x,b) = 
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For small values of b this diverges as log b for all values of x. Nonetheless, the SQM produces reasonable structure 
functions and quark distribution functions [20| . 

Another example is the NJL model, as regulated by two Pauli-Villars subtractions. The form factor is given by [n| 
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F^Q Z ) = / dxF T (Q*,x) (17) 
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with F^iQ 2 



— 12 ji ln(M 2 + A 2 + x(l — x)Q 2 ) TCg , where M is the quark mass, A is a parameter related 

to regularization and / is the pion decay constant. The subscript reg denotes the regularization procedure [22I ]: 
O rcg (A 2 ) = O(0) - C(A 2 ) + A 2 ^gr. The phenomenologically determined values are M = 280 MeV, /=93.3 MeV, 
A = 870 MeV. The impact parameter dependent pdf is the two-dimensional Fourier transform of F^(Q 2 ,x): 



q{x, b) 



-Ml 2 
(2tt)3/2 



dQQJ (Qb) 
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M 2 +x{\-x)Q 2 

M 2 + A 2 +x(l-x)Q 2 + M 2 + A 2 + x(l — x)Q 2 



This gives a well-behaved expression for b — > for all non-zero values of x(l 



Indeed: 
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Thus a logarithmic divergence appears upon integrating on x. 

Gaussian models with generalized parton distributions H(x,0,Q 2 ) (J dxH(x,0,Q 2 ) = F(Q 2 )) dominated by be- 
havior near x = 1 present a set of examples that also yield a form factor with a 1/Q 2 asymptotic behavior, and 
have a impact parameter distribution that is well behaved at each value of x for all b. The key asymptotic fea- 
tures are captured in the simple formula [H, HJ]: H(x, 0, Q 2 ) x ^i = (1 - x ) n ~ 1 e~ a ( 1 ~ x ^ Q , n > 2 so that 
q(x, b) x ^\ — 2 -na(i-x) e ~ b ^ iat " 1 ^ x ^ n ^ ■ This form shows that q(x,b) is well behaved for all values of b and for each 
value of x, but the integral over x contains a logarithmic divergence. 

Not all models that describe the exis ting form factor data have a singular central charge density. Relativistic light- 
front constituent quark models [25|, 12s [27j are able to describe the pion phenomenology and the current form factor 
data. These models produce a non-singular transverse charge density as we shall illustrate. These models can be most 
simply derived '26j by using the impulse approximation (evaluating the triangle diagram). One starts by evaluating 
the integral over the minus component of the loop momentum k^, and then cutting off the remaining integral over 
x = k + /p + , k±_ using a phenomenological wave function that depends on the combination (fc 2 + m 2 )/x(l — x), with 
m as the assumed constituent quark mass. We illustrate these models by computing the form factor using the model 
of [27] ]. The wave function chosen to be a power-law form, and the model is able to describe all of the existing form 
factor data in both the time- like and space like regions, f n , and the transition form factor f-^-y in which a virtual 
photon transforms a real pion into a real photon. The model form factor of [27J and the monopole fit of Eq. ([5]) are 
shown along with the measured data in Fig. 1 . Both models provide a good fit to the data, but present very different 
predictions for larger values of Q 2 where measurements remain to be done. The corresponding versions of p(b) and 
bp(b) are shown in Fig. 2. The singularity contained in Eq. @ appears as a rapidly rising function as b approaches 
zero, while the relativistic constituent quark model provides a p{b) that is smooth for small values of b. 

We summarize. The high Q 2 behavior of the form factor determines the short distance behavior of the transverse 
charge density p(b). If the form factor really behaves as the monopole form of E q . (H p) , then p(b) maintains a 
logarithmic singularity at the origin. A variety of models predict this behavior, [l3L Il4l. flEl Il9l I20L 123. I24]] as well as 
any non-relativistic constituent quark model that predicts a monopole behavior of F^iQ 2 ) and solutions of the Dirac 
equation. 

If the form factor falls asymptotically as perturbative QCD predicts, the p(b) behaves singularly as In In b for small 
values of b. It seems reasonable that p(b), a property of the valence quark density, should have no singularity. The 
relativistic constituent quark model produces transverse charge densities that are free of singularities, while providing 
a generally good phenomenology of the pion 27 1. 

It is therefore absolutely and manifestly clear that obtaining data at higher values of Q 2 is essential to providing 
further understanding. Such data could provide support for or rule out either constituent quark models or current 
pQCD evaluations of F^. If an assumption that the central density is non-singular is correct, the form factor will fall 
as described by constituent quark models. On the other hand, if asymptotic pQCD is valid, the central charge density 
would be singular- a remarkable fact of nature. 
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experimental data and the Physica file used for Fig. 1 . 
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